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We use the standard, adiabatic shell evolution to predict the size distribution N(R) for pop- 
ulations of SN-driven superbubbles in a uniform ISM. We derive N(R) for simple cases of 
ON . superbubble creation rate and mechanical luminosity function. We then compare our predic- 

tions for N(R) with the largely complete H I hole catalogue for the SMC, with a view toward 
the global structure of the ISM in that galaxy. We also present a preliminary derivation for 
N(v), the distribution of shell expansion velocities. 
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1. Introduction 

Core-collapse supernovae (SNe) tend to be correlated in both space and time because 
of the clustering of the massive (> 8M©) star progenitors. These clustered SNe, along 
, with stellar winds of the most massive stars, produce superbubble structures in both the 

warm ionized (10 4 K) and atomic H I components of the interstellar medium (ISM) in 
star-forming galaxies. The hot, coronal component of the ISM is thought to originate 
largely from the shock heating of material interior to shells of superbubbles and supernova 
remnants (SNRs). Total kinetic energies deposited into the interstellar environment are 
in the range 10 51 — 10 54 erg for OB associations, and > 10 55 erg for starburst phenomena. 
Hence, the large-scale structure and kinematics of the multi-phase ISM could be largely 
determined by this superbubble activity. Likewise, this effect should influence turbulence 
on global, macroscopic scales, which then cascades to smaller scales. 

The standard model for understanding superbubble evolution is the adiabatic model for 
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a continuous mechanical energy input (Pikel'ner 1968; Weaver etal. 1977; Dyson 1977), 
where the sequential SNe from the parent star cluster arc treated as a continuous power 
source (McCray & Kafatos 1987). The evolution of the shell parameters can then be 
described by a set of simple, self-similar relations analogous to the Sedov (1959) model 
for a single point energy injection. How applicable is the standard, adiabatic model 
to the long-term evolution of superbubbles? It is possible to make some rudimentary 
assumptions about the global ambient ISM and shell creation history, and then use the 
analytic equations to derive a superbubble size distribution that is predicted by the 
model. Such a size distribution can then be compared to observed H I shell catalogues 
to gain insight into the global structure and kinematics of the ISM. Variations from the 
prediction can then point to important effects that have not been adequately treated. 

In this contribution, we summarize our derivation of the superbubble size distribution, 
which is described in greater detail by Oey & Clarke (1997). We compare our results 
to the Hi shell catalogue of the Small Magellanic Cloud (SMC). Finally, we present 
preliminary new results for an analogous derivation of the distribution of shell expansion 
velocities. 



2. Assumptions 

Our purpose is to make the simplest feasible assumptions to see what the standard 
shell evolution predicts in the simplest conditions. Our assumptions are thus as follows. 
We assume coeval star formation in the parent clusters, which produce a constant 
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mechanical power L until the lowest-mass SN progenitors expire at an age t c = 40 
Myr. Studies of OB associations in the Magellanic Clouds and the Galaxy (Massey 
et al. 1995a, b) show stellar age spreads of < 3 Myr, motivating our adoption of coeval 
star formation. Likewise, stellar population synthesis modeling by Shull & Saken (1995) 
and Leitherer & Heckman (1995) suggests that the assumption of constant L appears to 
be reasonable. We also assume that the stellar initial mass function (IMF) remains fixed 
and universal, and we adopt a uniform and infinite ambient ISM. 

Based on the observed power-law form of the H n region luminosity function (H n LF; 
e.g. Kennicutt etal. 1989), we infer a similar power-law distribution of parent cluster 
masses. For a constant IMF, the nebular Ha luminosity scales directly with the number 
of stars, and likewise, L for the clusters is directly proportional to the number of core- 
collapse SN progenitors. We therefore take the mechanical luminosity function (MLF) 
for the cluster population to be a power-law as well: 

4>(L) = ^;=AL-e , (2.1) 

normalized such that J <fi(L) dL = 1. Ordinarily, the power-law index f3 of the MLF 
should be identical to that of the associated H II LF. We caution that evolutionary effects 
and small- number statistics in the stellar population can complicate this assumption (Oey 
& Clarke 1997, 1998), but essentially the power-laws are the same. We also consider one 
scenario with a single-valued MLF. In conjunction with these forms of the MLF, we 
consider a constant shell creation rate tp, and a single-burst creation model. 

The treatment of the endstage evolution for the superbubbles is crucial, but extremely 
uncertain. We assume that the shell growth stalls at an age t{ when the superbubble's 
internal pressure P\ < Pq, the ambient ISM pressure. Such a scenario is supported by 
numerical simulations (Garcfa-Segura & Franco 1996), in which radiative energy loss 
at this endstage suppresses further growth of the superbubble cavity. We then assume 
that the shell maintains this stall radius R{ until the input power stops at time t e . 
However, objects that never achieve pressure equilibrium with the ambient ISM continue 
to grow until t e . The subsequent destruction of the shells is even more uncertain. We 
simply assume that all objects survive for a nominal, universal period t s <C t c and vanish 
thereafter. We note that if the breakup of shells into smaller holes occurs such that 
the ratio of subunit sizes is universal for all objects, then an original power-law size 
distribution remains unaffected (Clarke 1996). 



3. Analytic Prediction 

The standard evolution predicts that the shell radius grows as (e.g. Weaver et al. 1977): 

/^50_\ 1/J 1/5 1/5 3/5 

V3087ry H V ' 

where t is elapsed time, p is the mass density of a uniform ambient medium, and all units 
are cgs. The pressure interior to the shell will decline as: 

p L^ L 2/5 3/5^4/5 (33) 

(3850tt) 2 / 5 f v ; 

From equations |3.2| and |3.3| we see that the stall criterion P; = Pq yields a correspon- 
dence between stall age and radius, uniquely determined by the input L. The charac- 
teristic time scale t e therefore determines the associated parameters R e and L c ; hence, 
a superbubble with input power L e will stall at exactly age t c with radius R c . It is thus 
apparent that objects with input power L < L c will follow an evolution that stalls at 



M.S. Oey & C.J. Clarke: The Superbubble Size Distribution 



3 



some point before t e , and remain at radii R < R e ; whereas those with L > L c will never 
stall, and at some point before t e will grow to radii R > R c . For an ambient number 
density n = 0.5 cm~ 3 , mean particle weight fj, = 1.25, and Pq = 3 x 10~ 12 dyne cm -2 , 
the adopted t c = 40 Myr implies R c = 1300 pc and L c = 2.2 x 10 39 ergs" 1 . 

These characteristic parameters are useful as scaling parameters, hence we have, 



R /L\ 1/5 /^ 3/5 



Rc V L c I V t, 



Pi / L\ 2/5 /t^- 4/5 



(3.4) 
(3.5) 



1/2 

(3.6) 



Pq V L c J \t c 

and 

t c _ R ( _ ( L 
t c R c \ L c 

We now derive the differential superbubble size distribution N(R) for specific combi- 
nations of shell creation history and MLF. We define N(R) dR as the number of objects 
with radii in the range R to R + dR. 

3.1. Continuous Creation, Single Luminosity 

For a continuous and constant superbubble creation rate 4> and a single- valued MLF with 
4>{L) — Lq, the size distribution for growing shells is given by, 

Ngiow ( R )=^t — \ . (3.7) 

The size distribution for the stalled objects is clearly a (^-function at the stall radius 
associated with L Q , whose magnitude is determined by the length of the creation period: 

N staU (R) = 1>(U- tf(Lo)) ■ 5(R - Rt(L )) , (3.8) 

where tf(L ) and Rf(L ) are the stall parameters for a shell powered by Lq. The dis- 
tribution in R is therefore determined exclusively by the growing objects. Applying the 



relations for the standard evolution given above, equation 3.7 gives, 
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We define the power-law slope of the size distribution as a, analogously to that of the 



MLF (equation 2.1), such that N(R) oc R~ a . The case here is the only one for which we 
derive a positive power-law index in R, yielding —a — The positive index is induced 
by the single- valued MLF. 

3.2. Single Burst Creation, Luminosity Spectrum 
We now consider the inverse case of instantaneous creation of all the objects, with a 



power-law MLF given by equation 2.1, Here, the size distributions for the growing and 
stalled objects are given by, 

dR\~ 1 



N gmw (R) = N h <f>{L) [ — ) (3.10) 



and 

N^„„(R) = Nu rh(T.) ( 

dL 



N ataU (R) = N h d>(L) 1^1 , (3.11) 
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respectively, where Nb is the number of objects created in the burst. Applying the 
standard evolution, we obtain, 

j 1-/3 f D \ 4-5/3 , . x -3+3/3 

N grow (R) = 5AN h (l-F st ) f r J (3-12) 



and 



r 1-/3 / n \ 1-2/3 

iVstaiiCH) = 2AN h F st — J , (3.13) 



where F s t is the fraction of stalled shells: 

Fst - ^ AZT' 3 <1L , (3.14) 

where L st (t) is the luminosity corresponding to t = if, i.e., the largest stalled shells. The 
lower limit of integration, £ m in, is the lower- L cutoff in the MLF, which in our analysis 
corresponds to the mechanical power associated with individual SNc. 

The HnLF in nearby galaxies typically has a power-law index of 2.0 ± 0.3 (e.g. Ken- 
nicutt etal. 1989), implying that value for j3. Hence, the power-law exponents of N(R) 
given by equations 3.12 and p. 13 are negative, and we have a ~ 6 and 3, for growing and 



stalled objects, respectively. 

We find that the stalled superbubbles generally dominate the total N(R), owing to 
the large numbers of weak-L objects. Therefore, the observed shell size distribution 



should be described essentially by equation 3.13, having a slope —a = 1 — 2/3. However, 



this will extend only to the stall radius associated with the age of the burst th, since 
higher-i objects will not have had enough time to stall. Objects with R > i?f (it>) must 
therefore be growing shells, described by equation 3.12| , having a much steeper dropoff 



with — a — 4 — 5/3. There is a discontinuous jump in N(R) by a factor of | at Rffa)- 

We therefore suggest that the existence, and age, of a recent single burst event can be 
discerned from the shell size distribution. 



3.3. Continuous Creation, Luminosity Spectrum 

In the most general case, continuous creation with a power-law MLF, the two different 
shell evolutions for L < L c and L > L c yield derivations for N(R) that produce different 
solutions for the regimes R < R c and R > R c . Since R c is very large (> 1 kpc), we are 
primarily interested in the regime for R < R e . 

For this case, superbubbles in the growing phase are described by, 

rU(R) /'dR\~ 1 
N gmw (R) = / V </>(£) 7J7 dt ■ ( 3 - 15 ) 

The lower limit of integration corresponds to the youngest objects having radius R, 
which are set by an upper-L limit L max . The upper limit of integration is the stall age 
corresponding to R (equation |3.6j ). For the standard evolution, equation 3.15 yields, 

7-1-/3 t / R\ 2 - 2P 
AT grow (i?) = BAi/j ^— — ^— — . (3.16) 



R c -2 + 3/3 \R, 
Similarly, stalled objects are given by 



X,„,i(/.'! = / IxftL) { dt . (3.17) 



dL 
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Applying the standard evolution: 
N staU (R) = 2Ai> 



1-2/3 



1 



Finally, objects surviving for a period t s , beyond t e 



{R) =L 



*c+t 3 



i> dt 4>(L) ( 



R 

R c 

are described by, 
dRf 



yielding, 



N sm {R) = 2Ayj 



L 



1-/3 



R c 



1-2/3 



(3.18) 



(3.19) 



(3.20) 



Adding together equations |3.16[ 3. IS, and 3.20| for the populations in the three evolu- 



tionary stages, we obtain an overall size distribution (for (3 > §): 



N(R) = Ai> 



LI"? 
R c 



1-2/3 



2(t c + t s ) + 



9-6/3 



U — 



(3.21) 



-2 + 3(3 

Equation p. 21 has two terms in R. As in the single-burst case, the stalled shells with the 
dependence —a = 1 — 2/3 dominate N(R). An observed shell size distribution would 
therefore resemble this power- law in the range 7?f(L m i n ) < R < R c . The lower limit 
represents the smallest stalled shells, which, in our treatment, correspond to individual 
SNRs. We caution that the assumption of constant L breaks down in this regime, since 
individual SNe are discrete events. However, the peak in N(R) should still be due to 
these single SNRs. 

For R < i?f(L m j n ), we recover growing objects, for which equation 3.15 is now domi- 
nated by the lower limit of integration, yielding: 

" -0 _ t / T? \ 2 /3 



N grov! {R) = 5A*P 



Li-? a 



Lc 



-2 + 3/3 



(3.22) 



Note that for this population, we again find the positive power-law slope —a = |, 
that we obtained for a single- valued <p(L) (equation |]^). The case here is analogously 
dominated by the upper-i limit, L max . 

Th e size distr ibutio n in the regime R > R c is derived from relations similar to equa- 
tions 3.15 , [3.17 , and 3. IS , with different limits of integration. We derive a final size 
distribution for these supergiant shells: 

Li-? 



N(R) = hAip 



R,. 



4-5/3 



-2 + 3/3 



(3.23) 



This population consists entirely of growing objects, along with a few shells in the survival 
stage. As seen in the case of the single burst, we again find a steep size distribution 
with —a = 4 — 5/3. Of course, given the numbers of objects, and breakdown of basic 
assumptions at these sizes, it is impracticable to compare the result to observations in 
this regime. 



4. Comparison with Observations 

Now taking the observed slope of the HuLF to be the same as /3, as argued in § 2, 
we can make a prediction for N(R) and compare with observed H I shell catalogues. 
Figure 1 shows the histogram for the size distribution of H I holes observed in the SMC 
(Staveley-Smith etal. 1997). We predict N(R) for this galaxy using the HllLF slope 
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Figure 1. Histogram of catalogued Hi shell radii in the SMC (Staveley-Smith etal. 1997). 
The observed power-law slope a Q = 2.7 ± 0.6 was fitted from the bins with error bars (dashed 
line); the predicted slope a p = 2.8 ± 0.4 is shown by the solid line, normalized to the data at 
R — lOOpc. The H I survey resolution limit is shown by the vertical long-dashed line. [Note 
that the slope of logiV(log R) vs. logi? is 1 — a.] 



measured by Kennicutt et al. (1989). We fitted a power-law slope to the H I data (dashed 
line) from the bins marked with error bars, yielding the observed slope a = 2.7 ± 0.6. 
The slope of our prediction (solid line) a p = 2.8 ± 0.4 is computed using (3 = 1.9 ± 0.2 
from the HilLF, and is therefore completely independent from a a . We normalized the 
prediction to the observations at R = 100 pc. 

The shell catalogue for the SMC is by far the most complete available for any galaxy. 
We also examined M31, M33, and Holmberg II (see Oey & Clarke 1997), but the H I data 
for those galaxies is too incomplete to discuss in detail here. For the SMC, however, the 
ratio of H I holes with R > 100 pc, to H n regions with Ha luminosities > 10 37 erg s -1 , is 
consistent with the life expectancies of the holes and nebulae, assuming constant creation. 
We are therefore confident that the SMC shell catalogue is reasonably complete at those 
radii. Compilation of an updated version of the catalogue is in progress, in which the 
very largest radial bins are slightly augmented (Stanimirovic etal. 1998). As seen in 
Figure 1, these new data may strengthen the derived power-law at large R. 

The observations and prediction are in remarkably good agreement, considering the 
many simplistic assumptions made in § 2. We therefore suggest that the SN-driven 
superbubble activity is indeed the primary agent responsible for the SMC shells, and 
that no other fundamental process is necessary to explain the H I superbubble structure 
in this galaxy. It is interesting to note that the —a = 1 — 2/3 power-law appears to 
be fairly robust, since it results from any evolution for which equation 3.6 applies. In 
particular, the rival momentum-conserving shell evolution also follows this relation, hence 
yielding the same power-law exponent for -/Vstaii(-R)- 

It will be interesting to examine N(R) for more galaxies (e.g. Thilker 1998), and espe- 
cially disks, where shell evolution should be affected by the non-uniform gas distribution. 
Surprisingly, the limited Hi hole samples for M31 and especially, M33, do not exhibit 
evidence that a power-law N(R) is truncated by blowout effects (Oey & Clarke 1997). 



M.S. Oey & C.J. Clarke: The Superbubble Size Distribution 



7 



SMC 




0.0 0.5 1.0 1.5 2.0 

log v(km/s) 



Figure 2. Histogram of shell expansion velocities for the SMC shell catalogue (Staveley-Smith 
etal. 1997). As in Figure 1, the power-law slope of —2.9 ± 1.4, fit from the data, is shown by 
the dashed line. The solid line shows the predicted slope of —3.5, normalized to the data at the 
last bin. 



5. Velocity Distribution 

Here, we briefly present a preliminary derivation of N(v), the distribution of super- 
bubble expansion velocities. We will carry out a more complete study in a future paper. 
We consider the case of constant -0 and power-law <j)(L), for the population of objects 
with R < R e . N(v) is clearly determined only by the growing objects, hence by analogy 
to equation |3.15 , 

rt(L max ,v) 

N gmw (v) = / V 4>{L) 

Jt(L min ,v) 



dv 
dL 



dt 



(5.24) 



where 



dR 
~dt 



3 R c 

5~r 



1/5 



-2/5 



(5.25) 



The limits of integration in equation 5.24 are the ages corresponding to the least and 
most luminous objects yielding v. We obtain: 

3/2-/3 



N(v) = -Aip 



25 



4-5/3 



5-5/3 



to 



-L 



1-/3 



.15/2-5/3 v - 



■7/2 



3-2/3 

(5.26) 

where v c is the soundspeed of the ambient ISM. This velocity distribution applies to 
objects that have not stalled, namely, those with v > v c . Th us, we find a power-law 
dependence of v~ 7 / 2 , independent of (3, although equation 5.26 is valid provided /3 > |. 

In Figure 2 we show the histogram of expansion velocities for the SMC H I shell cata- 
logue (Staveley-Smith etal. 1997). The fitted power-law slope is —2.9 ± 1.4, which is in 
agreement within the large error. The large uncertainty is caused primarily by the short 
dynamic range in logv, but Figure 2 shows that the power-law appears encouragingly 
well-behaved in comparison with the predicted slope. 
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6. Conclusion 

We have presented analytic derivations of characteristics of the superbubble population 
under various simple conditions, as predicted by the standard, adiabatic shell evolution. 
This approach is useful in determining the roles of various phenomena in the global 
structure of the ISM and evolution of superbubbles. The preliminary agreement found 
between our prediction and H I observations for the size and velocity distributions of shells 
in the SMC suggests that SN-driven superbubble activity is likely to be the dominant 
source of structure in the neutral ISM of this, and other, galaxies. The velocity structure 
and turbulence in the ISM are also likely to be substantially determined by this activity. 

Our analysis also yields other useful features that are discussed at greater length by 
Oey & Clarke (1997). For example, we predict a peak in N(R) at the stall radius of 
individual SNRs. Observations of N(R) might therefore provide an estimate for this 
radius, which could then be exploited to probe ISM densities, SN energies, and/or SNR 
evolution. The contribution of Type la SNRs should also be evident in N(R). Our 
analysis is also readily applicable to the porosity of the ISM, which determines the 
relative importance, by volume, of the hot, coronal gas compared to the cooler phases of 
the ISM. Our study derives analytic expressions for both 2D and 3D porosity parameters. 
For three of the four galaxies examined, we found porosity parameters < 0.3, suggesting 
a fairly low filling factor for the hot ISM, although we caution that these values are 
sensitive to assumptions for the ambient interstellar conditions. We also applied this 
porosity analysis to the Galaxy, with inconclusive results. 

We are grateful for discussions with E. Blackman, D. Cox, L. Drissen, J. Franco, D. 
Hatzidimitriou, C. Leitherer, C. Robert, J. Scalo, J.M. Shull, L. Staveley-Smith, and G. 
Tenorio-Tagle. 
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